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The evolution properties of propagating particles produced at high 
energies in a randomly distributed environment are studied. The fi- 
nite size of the phase space of the multiparticle production region as 



c~| , well as the chaoticity can be derived. 
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1 Introduction 

Interest in (charged) particles "moving" in an environment of quantum fields 
taking into account the relations between quantum fluctuations and chaotic- 
ity is still attracted by the particle physics society. The events of high multi- 
plicity at high energies is an interesting subject because of increasing energies 
of current and future colliders (e + e~, pp, pp) of particles. One of the most 
important tasks of superhigh energy particle studies is to analyze fluctuations 
and correlations such as the Bose-Einstein (BE) correlation [1] of produced 
particles. This is a rather instructive tool to study high multiplicity hadron 
processes in detail. We understand the multiparticle (MP) production as the 
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process of colliding particles where the kinetic energy is dissipated into the 
mass of produced particles [2]. We consider the incident energy yfs ^> A, 
where A means the quantum chromodynamics (QCD) scale. Phenomeno- 
logical models [3,4] describing the crucial properties of multiparticle corre- 
lations are very useful for systematic investigations of the properties caused 
by fluctuations and correlations. By considering them, one can obtain the 
characteristic properties of the internal structure of disordering of produced 
particles in order to extract the information on the space-time size of the 
multiparticle production region (MPPR), to estimate the lifetime of the par- 
ticle emitter, etc. There were used the analysis of correlation functions (CF) 
and distriburion functions (DF) in [5,6] to understand the possible view on 
the quark-gluon plasma formation. In this paper, we present the model to 
describe the very high multiplicity effects at high energies. The most char- 
acteristic point of our model is that the BE both DF's and CF's are taken 
into account on the quantum level (the operators of production and anni- 
hilations are used) with the random source contributions coming from the 
environment. We define the average mean multiplicity (N) via MP CF w(k) 
as (N) = J d$kw(k), where k is the spatial momentum of a particle. Fol- 
lowing a natural way we suppose (N) <C N, while Af <C iVo = y/s/m, where 
m ~ 0(0.1 GeV). The main object in this investigation is the MP thermal 
DF W(k^ k^) related to (N) as 

W{k») = (N)f{kJ = (N)(b + (k,)b(k^) p , (1.1) 

where (N) is defined as the scale of the produced-particle multiplicity N 
at four-momentum k^ (/j, is the Lorenz index), the normalized DF f{k p ) is 
finite, i.e. / d^k f(k^) < oo and the label f3 in (|1.1|) ) means the temperature 
T (of the phase space occupied by operators b + (k fl ) and b{k^)) inverse. The 
nature of operators b + (k fl ) and b{k^) will be clarified in Sec. 3. 

In this paper, we claim that the observation of the size effect in a MP 
production is derived via MP CF's and DF's as well as the so-called chaotic- 
ity which will be introduced in Sec. 4. The MP CF formalism concerns the 
statistical physics based on the Langevin-type of equations. The Langevin 
equation to be introduced in Sec. 3 is considered as a basis for studying the 
approach to equilibrium of the particle(s). It is assumed that the heat bath 
being in essence infinite in size remains for all times in equilibrium as well. 
We use the method applied to the model where a relativistic particle mov- 
ing in the Fock space is described by the number representation underlying 



2 



the second quantization formulation of the canonical field theory. We deal 
with the microscopic look at the problem with the elements of quantum field 
theory (QFT) at the stochastic level with the semiphenomenological noise 
embedded into the evolution dissipative equation of motion. 

2 The model of dual representation within 
the dissipative dynamics 

Let us suppose that the evolution of particles produced at high energies is 
described by the solutions of the model Hamiltonian where any physical sys- 
tem of particles is described by the doublet of field operators. We introduce 
the dual states in the eigenbasis {\k), \k(e))} of the Hamiltonian H, where 
< e < oo is identified as the frequency representing the energy of the object 
and \k) is a discrete eigenstate. We consider the simple dual model where H 
is given within the damped harmonic oscillator 

H = H + H l} (2.1) 

where 

POO 

H = e \k)(k\+ / dee\k(e))(k(e)\ , (2.2) 
J o 



H t = p deg(e)[\k(e))(k\ + \k)(k(e)\], 
Jo 



(2.3) 



(k\k) = 1 , (k\k(e)) = (k(e)\k) = , (e|e'> = 5(e - e') . 

Here, we identify \k), \k(e)) and (k\, (k(e)\ with the special mode operators 
of annihilation a, b(e) and creation a + ,b + (e), respectively, e.g., for "quarks" 
and "gluons" or their combinations. In the interaction Hamiltonian ( |2.3| ), p 
is the coupling constant, while g(e) provides the transition between discrete 
and continuous states. 

The equations of motion obeying (|2.1| ) with (|2.2|) and fl2.3|) are 



i d t a k (t) = e a k (t) + p / de g(e) b k (e,t) 

Jo 

i d t b k (e, t) = e b k (e, t) + p g(e) a k (t) , 
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where the label k means \k\ as the momentum. We demand that ak(t) and 
b k (e, t) satisfy the following natural conditions: 

/■oo 

a k (t)a£(t)+ deb k (e,t) b~£ (e,t) < oo , 



[a+(t)]+ = a k (t) ,[b+(e,t)]+ = b k (e,t). 

The next step is to give the relations between the variables and couplings of 
the model (|2.1|) and the phenomenological constants, e.g., the frequency E 
and the damping constant k involved into the dissipative dynamics given by 
the Boltzmann-type equations in the relaxation time approximation: 

= -«* + «), (2.4) 
dtw k (t) = d t (at(t)a k (t)) = 2 k [w k (t) - n(f3)] . (2.5) 

Here, 

poo pt 

d t {a k (t)) = -ie (a k (t)) - p 2 deg 2 (e) J ds (a k (s)) exp[-ie(t - s)] , 
(a fe (z)) = -ip J dx (cr k (x)) exp[-ie (^ - x)} , 

poo 

&k{t) = / deg(e)b k (e,t) , 



o 



n(/3) = [exp(e /3)±l] 



-i 



as r — >■ — oo. One can conclude that the phenomenological constants E and 
k in ( |2.4| ) ( p.5|) are nothing else but e and 2 7rp 2 <? 2 (e ), respectively, i.e. the 



microscopic parameters coming from the model Hamiltonian (EO 



3 Stochastic model. Langevin equation 

As it was pointed out in the Introduction, to derive the characteristic fea- 
tures of the MP production physics at high energies, one should specify the 
model on the quantal level. Let us assume that only the particles pi of the 
same kind of statistics labeled by index i are produced just after the high 
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energy collision process occured, e.g., pp,pp — > pi. In order to extend the 
method of stochastically distributed particles in the environment, we pro- 
pose that rather complicated real physical processes happened in the MP 
formation region should be replaced by a single-constituent propagation of 
particles provided by a special kernel operator (in the stochastically evolu- 
tion equation) considered as an input of the model and disturbed by the 
random force F [5,6]. We assume that F can be the external source being 
both a c-number function and an operator. In such a hypothetical system of 
excited (thermal) local phase we deal with the canonical operator a(k, t) and 

— # 

its Hermitian conjugate a + (k,t). We formulate DF of produced particles in 
terms of a point-to-point equal-time temperature-dependent thermal CF of 
two operators 

w(k, k', t; T) = (a + {k, t) a(k', t)) p = Tr[a + {k, t) a{k', t)e- Hf} \/Tr{e- Hf3 ).{Z.l) 

Here, (...)/3 means the procedure of thermal statistical averaging; k and t are, 
respectively, momentum and time variables, e~ H " /Tr(e~ HI3 ) stands for the 
standard density operator in equilibrium, and the Hamiltonian H is given by 
the squared form of the annihilation a p and creation operators for Bose- 
and Fermi- particles, H = Y^p^Ap^p (the energy e p and operators a p , a+ 
carry some index p, where p a = 2tt n a /L,n a = 0, ±1, ±2, V = L 3 is the 
volume of the system considered). Here, we use the canonical formalism in a 
stationary state in the thermal equilibrium (SSTE), and a closed structural 
resemblance between the SSTE and the standard QFT is revealed. We define 
the thermal boson field as 

$bM = 4? \<i>M + , (3.2) 



where 



V2 

00^) = J d 3 kv k a(k,t) ,v k = 
<t> + {x») = I d 3 kv£a + (k,t) 



[(2tt) 3 A(A;)] 1 /2 
~ [{2n) 3 ^{k)} 1 / 2 



and A(k) is an element of the invariant phase volume. 
The standard canonical commutation relation (CCR) 



a(k,t),a + (k',t)\ ± = 5 3 (k-k') (3.3) 
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at every time t is used as usual for Bose (-) and Fermi (+)-operators. 

The probability to find the particles in MPPR with momenta k and kl in 
the same event at the time t is: 

R(k, k', t) = W(k, k', t)/[W(k, t) ■ W(k\ t)] , (3.4) 

where the MP DF W(k,t) in the simple version fluctuates only its normal- 
ization, e.g., the mean multiplicity (N) : 

Here, the one-particle thermal DF is defined as 

W(k,t) = {N). f(k,t) , 

defining the single spectrum, while 

W{k,k',t) = {N{N' - 5^)) ■ f{k,k',t) 

for i-and j-types of particles. Here, Sij = 1 if i = j and otherwise. DF's 
f(k,t) and f(k,k',t) look like (hereafter the label f3 will be omitted in the 
sense of Q and Q)) 

f(k,t) = (b+(k 1 t)b(k,t)}, 
f(k,k',t) = (b + {k,t) b + (k',t) b(k,t) b(k\t)) , 

where 

b(k,t) = a(k,t) + R(k,t) 

under the assumption of the random source-function R(k, t) being an oper- 
ator, in general. One can rewrite (|3.4j) in the following form 



R(k,k',t) = m /a '- Kn 



/(M)/(^)' 

where 

urn = 

qy } (N)(N') ' 

For simplicity, we deal with operators a and b as if they are the single boson 
or fermion operators. Considering the "propagation" of a particle with the 
momentum k in the quantum equilibrium phase space under the infuence of 
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a random force coming from surrounding particles, the dissipative dynamics 
of the relevant system is described by the equation containing only the first 
order time derivatives of the dynamic degrees of freedom, the operators b(k, t) 
and b+(k,t) [5]: 

i d t b(k,t) = F(k,t) - A(k,t) + P , (3.5) 

% d t b + {k, t) = A*{k, t) - F + {k, t)-P . (3.6) 

Here, the interaction of particles under consideration with surroundings as 
well as providing the propagation is given by the operator A(k, t) defined as 
the one closely related to the dissipation force: 

A(k,t)= K(k,t-r) b(k,r) dr . (3.7) 

J — oo 

The particle transitions are provided by the random source operator F(k, t) 
while P stands for a stationary external force. An interplay of particles with 
surroundings is embedded into the interaction complex kernel K(k,t), while 
the real physical transitions are provided by the random source operator 
F(k, t) with the zeroth value of the statistical average , (F) = 0. The random 
evolution field operator K(k, t) in ( p.7|) stands for the random noise and it is 
assumed to vary stochastically with a 5- like equal time correlation function 

(K + (k, t) K(k\ t)) = 2 (7rp) 1/2 k 5{k - k') , 

where both the strength of the noise k and the positive constant p — > oo 
define the effect of the Gaussian noise on the evolution of particles in the 
thermalized environment. 

The formal solutions of (|3T5|) and ( |3.6|) in the operator form in S , (3? 4 ) 
(Jfef = (uj = k°,kj)) are 

6(/c M ) = a(*g + R(k^) , 

b + (k,)=ti + (k,) + R*(k,) , 

respectively, where the operator a(k^) is expressed via the Fourier trans- 
formed operator F(k^) and the Fourier transformed kernel function K(kfj) 
(coming from (|3.7|)) as 

a(k ll ) = F(k fl )-[K(k fl )-u}- 1 , 
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while the function R(k^) ~ P ■ [K(k^) — oj]^ 1 . In our model, we suppose that 
a heat bath (an environment) is an assembly of damped oscillators coupled 
to the produced particles which in turn are distributed by the random force 
F(k^). In addition, there is the assumption that the heat bath is statistically 
distributed. The random force operator F(k,t) can be expanded by using 
the Fourier integral 



F(k,t) 



(3.8) 



where the form ^{k^-cik^ is just the Fourier operator F(k^) = ip(k^)-c(k 
and the canonical operator c(k^) obeys the commutation relation 



HJi 



c(^),c + (fc' 



± 



$ 4 (k, - k') 



The function ipfe^) in (|3.8| ) is determined by the condition (the CCR ( |3.3| ) 
is taken into account) 



r+°° duo 


[ 1>%) 


Loo 2tt 





n 2 



4 DF ratio enhancement. 



The enhanced probability for emission of identical particles is given by the 
ratio R of DF's in Sffiz) as follows: 



R(k„ k' ; T) = £(N) 4 



fv^n-i hp] T) 



f(K) ■ f(K) 



(4.1) 



where /(*„, T) = (&+(*„) and f(k,) = (&+(*„) &(*„)). 

Using Fourier solutions of equations ( |3.5| ) and ( |3.6| ) in 5 , (3 f ?4) , one can get the 
R-ratio for DF 



R(k„ K; T) = £(iV) [1 + D{k„ k'; T)] 



(4.2) 



where 



n/ , y m\ _ Z{k^)[~{k>^) + R + {k'^)R{k,)] + ~(k^k,)R+(k,)R(k' fl \ A 

/(U • f{k' ) 
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and the Bose-Einstein CF k' ) looks like 

- ■ (c + (A; M ) c(Ag) . (4.4) 



[K*{k,)-u]-[K{k'\-u'} 



Inserting CF ( |4.4| ) into (|4.3|) and taking into account the trick with S 4 (k fl — 
A^J-function to be replaced by the 5-like consequence like fl(r) ■ exp[— (k — 
k') 2 r 2 } [7], one can get the following expression for D-function instead of (|4.3| ) 



D(k„ fcj,; T) = \(k„ T) • exp(-g 2 /2) 
x [71(0;, T)fi(r) exp(-g 2 /2) + R*(k'jR(k„) + R*(k„)R(k'J] , (4.5) 

where 

X(k u , k'\ T) = — — = ■ n(u, T) , d) = -iu + u') , 



while g 2 = Q 2 r 2 and the function Q(r)-n(u; T) -exp(— q 2 /2) in ( |4.5| ) describes 
the space-time size of the multiparticle production region. Choosing the z- 
axis along the pp or pp collision axis one can put 

Qn = (k — k )n, Qo = tj: — eg/, Q z = k z — k z , Q t = [(k x — k x ) + (k y — k y ) ] , 

^( r ) = r o • r z ■ r 2 , 

where ro, r z and r t are time-like, longitudinal, and transverse "size" compo- 
nents of MPPR. To derive (|4.5[) , the Kubo-Martin-Schwinger condition ( fi 
is the chemical potential) 

(a(k\ f) a + (k, t)) = (a + (k, t) a(k\ t - i(3)) ■ exp(-/5 fi) 

has been used, and the thermal statistical averages for the c(/c M )-operator 
should be represented in the following form: 

(c + (k,)c(k' ll )) = 5\k ll -k^-n(u J ,T) , (4.6) 
<c« c + {k' )> = 8\k, - k' ) • [1 ± n{u>, T)] (4.7) 
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for Bose (+)- and Fermi (-)-statistics, where n(u, T) = {exp[(u; — /i)/3]±l} _1 . 
Formula (|4.5| ) indicates that the chaotic multiparticle source emanating from 
the thermalized MPPR exists. Taking into account Q4.6|) and ( |4.7| ), it is 
easy to see that the correlation functions containing the random force func- 
tions F(k,t) (|3.8| ) carry the quantum features in the termalized stationary 
equilibrium, namely 



r^-At) = I —\ij(k^\ 2 [l±n(u,P)} exp(-iujAt) 



{F(k, t) F + (k', t')) = 5%k - k')) T^k, -At) 
duj 

At = t-t'; 
(F + (k, t) F + (k', t')} = 5 3 (k- k')) T 2 (k, At) , 

T 2 (k,At) = I ^\ip(k J\ 2 n(u,P) exp(iu;At). 
Jin 

The quantitative information (longitudinal r z and transverse r t components 
of MPPR, the temperature T of the environment) could be extracted by fit- 
ting the theoretical formula (3~5) to the measured D-function and estimating 



the errors of the fit parameters. Hence, the measurement of the space-time 
evolution of the multiparticle source would provide information of the mul- 
tiparticle process and the general reaction mechanism. The temperature 
of the environment enters into formula (|4.5|) through the two-particle CF 



H(fc M , k' ; T). If T is unstable, the /^-functions ( |4. 1| ) will change due to a 
change of DF / which, in fact, can be considered as an effective density of 
the multiparticle source. Formula (|4.1|) looks like the fitting R-ratio using a 
source parametrization: 

R F (r) = const [1 + X F (r) ■ exp(-rf ■ Q 2 t /2 - r\ ■ Q 2 z /2)} , 

where r t (r z ) is the transverse (longitudinal) radius parameter of the source 
with respect to the beam axis, stands for the effective intercept parameter 
(chaoticity parameter) which has a general dependence of the mean momen- 
tum of the observed particle pair. Here, the dependence on the source lifetime 
is omitted. The chaocity parameter is the temperature-dependent and 
the positive one defined by 

X F (r) W^^l 2 



f(K) ■ f(K) 
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Comparing ( |4.3| ) and ( [4.4| ) one can identify 



E(k„ fcp = O(r) • n(cu; T) • exp(-g 2 /2) . 

Hence, CF fc' ) defines uniquely the size r of MPPR. There is no satisfac- 
tory tool to derive the precise analytic form of the random source function 
R(k^) in ([Of), but one can put (see (|4.4j ) and taking into account R(k^) 
~ P ■ [K(k^ - to]- 1 ) [8,5,6] that 

R{k,) = [ a -~{k^)] l l\ 

where a is of the order O [P 2 j 'n(uj , T) ■ {^(k^f^j . Thus, 

D(q 2 ; T) = jS^T-I-^ 12 P V2 (^ T ) e " ?2/2 + 2 (^') 1/2 ] , (4-8) 
(1 + a){l + a') 1 J 

where 

A(u;T) 



n(uj; T) ■ n(u/; T) 

It is easy to see that, in the vicinity of q 2 ~ 0, one can get the full correlation 
if a = a' = and X(uj;T)=1. Putting a = a' in ( |4.8|) , we find the formal 
lower bound on the space-time dimensionless size of MPPR of the bosons 

2^1 A(o);T) 
q > in 



[^(a + l) 2 + a 2 - a] 



In fact, the function D(k fl , k';T) in ( |4.5|) could not be observed because of 
some model uncertainties. In the real world, D-function has to contain back- 
ground contributions which have not been included in the calculation. To be 
close to the experimental data, one has to expand the D- function as pro- 
jected on some well-defined function (in S^R.^)) of the relative momentum of 
bosons produced D(k^, k'^T) — > D(Q 2 l ;T). Thus, it will be very instructive 
to use the polynomial expansion which is suitable to avoid any uncertainties 
as well as to characterize the degree of deviation from the Gaussian distribu- 
tion, for example. In (— oo, +oo), a complete orthogonal set of functions can 
be obtained with the help of the Hermite polynomials in the Hilbert space 
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of the square integrable functions with the measure dfi(z) = exp(—z 2 /2)dz. 
The function D corresponds to this class if 

/+oo 
dqexp{-q 2 /2) \D{q)\ n < oo , n = 0, 1, 2, ... . 
-00 

The expansion in terms of the Hermite polynomials H n (q) 

D(q) = A £ c n - H n (q) ■ exp(-q 2 /2) (4.9) 

n 

is well suited for the study of possible deviations both from the experimental 
shape and from the exact theoretical form of the function D ( |4.5| ). The coef- 
ficients c n in (}4.9| ) are defined via the integrals over the expanded functions 
D because of the orthogonality condition 

/+oo 
H n (x) H m (x) exp(-x 1 2) dx = 8 n>m . 
-oo 

Thus, the observation of the multiparticle correlation enables one to extract 
the properties of the structure of q 2 , i.e. the space-time size of MMPR. The 
other possibility is related with the replacement of R- function ( [4.2| ) with 
respect to the cylindrical symmetry angles 6 and <fi which are non-observable 
ones at fixed Q t : 

R(k^ fc;; T) - R(Qu T) = C N X £(N) J dq t dQ z dd d<P f{k^ T) , 
where 

C N = I dq t dQ z dO d<f) f(k^) f(k'^) , 

Qt = CQS 6> + sin 6> { kx + ky T 2 Qt [C0S( ^ + 4>) + Sin( ^ + 0)] } • 
Then, R{Q t - T) = £(N) [1 + D(Q t ; T)) with 

5(ftir)= (i^Ki^-) exp| - (r?Q?)|f(ft;r) ' 

F{Q t ;T) = J dq t dQ z d8d<f)n 2 {u;T) e _fe [1 + 2 ^ aa'\-\u- T) e q2/2 } , 

p 0z = r 2 Q 2 + r 2 z Q 2 z , 
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Cjv(T) = dq t dQ z d9 d<f>n(uj;T) n(uj';T) . 



It remains to sketch how one goes about calculating the thermodynam- 
ical quantities in a local thermalized system of produced particles. Tak- 
ing into account the positive- and negative-frequency parts of the boson 
field operator ( |3.2| ) to be applied to the energy-momentum tensor T^ u (x) =: 

®b(x) k^k u : and the particle flow operator Il M (x) =: &b{x) k^ &b(x) '■ 

we can calculate the energy density E(/3), the pressure V(/3) and the entropy 
density S(/3) in the local system of the volume v for bosons in the equilibrium 
thermalized phase space. The simple straightforward calculations give (see 
also [9]) 

E ^) = TTW~ / d 3 kdcucu 2 M(k„P), 

(2 TV) V J 



D 7T V 



S(J3) = I d 3 k{[l + U (k,(3)] ln[l + U (k,P)]-U (k,(3) lnU (k,(3)} 

1 n v .1 



where 



U (k,(3) = - / d 2 uM{k^(3), 

L J — oo 

d?.k = 



;2tt) 3 A(A;) 



=. 1 

$ B (x) k^ $ B (x) = -[$ B (x) {k^ B {x)) - {k^ B {x)) $ B (V>] . 

Here, we suppose that the thermalized MPPR is isotropic, and one can use 
the space-averaged operators normalized to the volume v, taking ensemble 
averages Q4.6Q and ( [Op . It is easy to see that both E{(3) and V(j3) tend to 
their maximum values with rising T, while the entropy S(f3) changes not so 
much essentially even if T — > oo. 
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5 Statistical distributions 



From a widely accepted point of view, at high energies, there are two chan- 
nels, at least, for multiparticle production where produced particles occupy 
the MMPR consisting of % elementary cells. These main channels 
direct channel assuming that all particles pj are produced directly within 
the quark (g)-antiquark (q) annihilation or the gauge-boson fusion, e.g., 
qq — > PjPj...] b) an indirect channel which means that the particles are pro- 
duced via the decays of intermediate vector bosons x* i n both heavy and 
light sectors in the kinematically allowed region, e.g., qq — > x*X* ■■■ PjPj — 
All the particles produced are classified by the like-sign constituents that are 
labeled as p + ,p~ , p°-subsystems, where p : \i , n , K... The mean multiplicity 
(N) and the mean energy (E) of the pj subsystem are defined as [3] 

W = EE rn 3 <f '> , 

j mj 

j m 3 

where €j is the energy of a p-particle in the jth elementary cell and Cj™^ 
stands for the probability of finding rrij p particles in the jth cell and is 
normalized as 

oo 

E t 3) = 1 • 

TOj=0 

In the direct channel, for charged produced mesons (N) is defined uniquely 
for a given (5 as 

(iV) = 2E[exp(e^)-l]- 1 , 

3 

while (E) is 

(E> = ^ = ^expM)-r 

Going into y-rapidity space in the longitudinal phase space with a lot of cells 
of equal size Sy, the energy €j should be expressed in terms of the transverse 

mass m t = \J (k t ) 2 + ({k t ) and m p are the transverse average momentum 
and the mass of a p-particle): 
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9j (s) = (VS + y/g=l) exp[-(j - 1/2) 5y\ . 
Here, the four-momentum of £>-particle is given as 



k^ = (y(kt} + mp cosh y, k t cos <p, k t sin ip, y (k t ) +m^smhy) } 

where the azimuthal angle of k t is in the range < <p < 2n. Our model 
produces an enhancement of R(Q,/3) in the small enough region of Q where 
R is defined only by the model parameter a and the mean multiplicity (N(s)) 
at fixed value of (3, namely: 



R{Q,P)~£{{N)){1 



[1 + aY 



X(Q, (3) + 2 a - (yJ\(u, (3) + a)Q 2 r 2 



It is clear that R(Q, /3)-function at Q 2 = 

R(Q,{3)~t{(N(s)})- 2- 



a 



1 + a 



cannot exceed 2 because of a ^ and £(iV(s)) < 1 even at large multiplicity. 
The Boltzmann behaviour should be realized in the case when a — > oo, i.e. 
the main contribution to the fluctuating behaviour of the R(Q, (3) -function 
should come from the random source contribution (see ( |4.5| ) and (|4.8D ). We 
found that the enhancement of the R(Q, (3) -function, mainly, the shape of 
this function, strongly depends on the transverse size r t of the phase space 
and has a very weak dependence of the Sy size of a separate elementary cell. 
The increasing of r t makes that the shape of the R(Q, /3)-function becomes 
more crucial. 

Obviously, £((iV(s))) is the normalization constant in (|4.2| ), where (N(s)) 
should be derived at the origion of Q 2 precisely from R(Q = 0, (3) = Rq(s) 
as 



where 



(N(s)) ~ - 



Ro(s) 



l+Q 



can be extracted from the experiment at some chosen value of a (a = 
should be taken into account as well). On the other hand, the R(Q,(3)- 
function allows one to measure a = a' which parametrizes the random source 
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contribution as well as the splitting between a and a'. Neglecting the random 
source contribution (i.e., putting a = a' = 0) we can estimate the chaoticity 
A(u>,/3) by measuring R(Q,{3) as Q 2 — > 0. 

In fact, the theoretical prediction that D(Q, f3) > 1 means that in MPPR 
one should select the single boson "dressing" of some quantum numbers, and 
the particles suited near it in the phase space are "dressed" with the same 
set of quantum numbers. The amount of such neighbour particles has to be 
as many as possible. This allows to form a cell in the space-time occupied 
by equal-statistics particles only. Such a procedure can be repeated while 
all the particles will occur in the MPPR. This leads to the space-time BE 
distribution of produced particles in the phase-space cells formed only for 
bosons. In fact, there is no restriction of the number of bosons occupying 
the chosen elementary cells. It means that D(Q, /3)-function are defined for 
all orders. 



6 Summary and discussion 

We investigated the finite temperature BE correlations of identical parti- 
cles in the multiparticle production using the solutions of the operator field 
Langevin-type equation in SiJR,^), the quantum version of the Nyquist the- 
orem and the quantum statistical methods. The model considered states 
that all the particles are produced directly from the a high-energy collision 
process. There was presented the crucial role of the model in describing the 
BE correlations via calculations of DF's as functions of the mean multiplicity 
and chaoticity at each four-momentum \JQ\\- Based on this model, one can 
compare the effects on single particle spectra and multiparticle distribution 
caused by multiparticle correlations. There are several parameters in the 
model: /3 , 8y , a(a') . One can focus on the statement that the deviation of 
the D(Q, (3) -function from zero at finite values of the physical variables q 2 
and the model parameter a indicates that the MPPR should be considered as 
the phase-space consisting of the elementary cells (with the size 5y of each) 
which are occupied by the particles of identical statistics. The Boltzmann 
behaviour of i?-function is available only at large enough values of a which 
means the leading role of the random source contribution to the distribution 
function. An important feature of the model is getting the information on 
the space-time structure of the multiparticle production region. We are able 
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to predict the source size and the intercept parameter -the chaoticity A as 
well. We have found that DF R(Q, j3) depends on the number of elementary 
defined by the size Sy in the rapidity ?/-space. 

Of course, the best check of any model could be done if various kinds of 
high energy experimental data on the multiparticle correlations would be well 
reproduced by the model in consideration. Our model can be applied to the 
data which are available from the ALEPH experiment [10] at \/s=91.2 GeV. 
The ALEPH data applied to _R(Q)-function (in the notation of our model) in 
the region O.lGeV < Q < lGeV have an essential nonmonotonic behaviour. 
It means that the sign of the slope parameter of R(Q) changes in this region, 
and the values of R(Q) are less than unity. In accordance with our model, 
this effect could be interpreted as the fact that in the domain O.lGeV < 
Q < lGeV the mean multiplicity (N(s)) has a small value which gives the 
strong suppression for the i?(Q)-function normalized to £((iV(s))). In the 
case of higher energies of colliding particles (Tevatron, LHC) , the expected 
(from our model) averaged multiplicity should be very large ~ 300-400, and 
the nonmonotonic effect will not occur. We hope that experiments at the 
Tevatron and LHC will measure the observables like R(Q) more precisely 
and the puzzle regarding the nonmonotonic behavior of the multiparticle 
distribution will be clarified. 
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